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We introduce a modified factorization formalism in quantum chromodynamics for hadronic 
production of W and Z bosons at large transverse momentum pr- When pt is much larger than 
the invariant mass Q of the vector boson, this new factorization formalism systematically resums 
the large fragmentation logarithms, ln’”(py/Q^), to all orders in the strong coupling as- Using 
our modified factorization formalism, we present the next-to-leading order (NLO) predictions for 
W and Z boson production at high pt at the CERN Large Hadron Collider and at a future 100 
TeV proton-proton collider. Our NLO results are about 5% larger in normalization, and they 
show improved convergence and moderate reduction of the scale variation compared to the NLO 
predictions derived in a conventional fixed-order perturbative expansion. 
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I. INTRODUCTION 

Successful operation of the CERN Large Hadron Col¬ 
lider (LHC) and associated particle detectors have led to 
the discovery of the Hi^s boson, the final piece of the 
standard model (SM) of particle physics, along with 
the exploration of strong interaction short-distance phe¬ 
nomena in high energy processes at much greater values 
of the production transverse momentum. Future experi¬ 
mental investigations, at higher collision energy and with 
greater luminosity, promise refined understanding of the 
nature of electroweak symmetry breaking and possible 
evidence of new physics beyond the SM. Some of these 
searches will focus on deviations from SM expectations 
or on anomalies in high-energy tails of various kinematic 
distributions. Precise SM predictions for these observ¬ 
ables and distributions are important assets for discovery 
of new physics. 

The production distributions of massive electroweak 
(EW) gauge bosons, W’s and Z’s, are among observables 
that can be sensitive to physics beyond the SM, either be¬ 
cause new states in extensions of the SM may decay into 
W and Z bosons, or because the lepton distributions from 
SM W and Z decay are important backgrounds for high- 
energy lepton signatures in new physics models. The lep¬ 
tons could mimic the signature of boosted objects from 
the decay of a new heavy resonance. Moreover, W and Z 
boson production serve as tests of perturbative quantum 
chromodynamics (QCD) calculations, and data on their 
distributions are important in the determination of the 
parton distribution functions (PDFs) Q. Among mea¬ 
surements of W and Z production, precise cross sections 
at large transverse momentum hold particular interest. 
The SM Z boson production at large pr can be used for 
jet energy-scale calibration. Existing studies show the 
possibility of using W and Z boson production at large 
Pt to further constrain the gluon PDFs [1, [^. Trans¬ 


verse momentum distributions of W and Z bosons have 
been measured at the LHC [1,[3| but with relatively large 
uncertainties. 

To exploit the full potential of the data, ever more 
precise SM predictions are needed, requiring better un¬ 
derstanding of the size of high-order perturbative correc¬ 
tions. The Pt spectra of W and Z boson production in 
hadronic collisions have been calculated perturbatively in 
the SM. The leading order (LO) cross section at a finite 
Pt is of 0{aemC(s)- The next-to-leading order (NLO) 
QCD corrections were calculated decades ago [a-[lOl|. and 
found to be sizable for the energy of the LHC. The EW 
corrections were studied extensively in recent years [n^ 
fisjl . The W and Z production at small pt Q mass 
of W or Z) has received a lot of attention in connec¬ 
tion with the resummation of the effects Sudakov double 
logarithms in QCD [T^-[3^, also important for precise 
measurements of the W boson mass in hadronic colli¬ 
sions. Improved predictions are awaited of W and Z 
production at large pt beyond NLO corrections [33 - [4ll| . 
Work is in progress on the full next-to-next-to-leading 
order (NNLO) QCD corrections to the pt spectrum of 
W and Z production, similar to the case of Higgs boson 
production 

In this paper, we explore another approach to theoret¬ 
ical improvement by identifying potentially large QCD 
logarithms from high order perturbative calculations, 
and resumming these logarithms to all orders. QCD cor¬ 
rections to the short-distance partonic scattering cross 
sections of W and Z production at large pt could receive 
one power of a large logarithm, hi{p^/Q^), for every ad¬ 
ditional power of Us- Such large logarithms come from 


^ After our paper was submitted, the NNLO QCD corrections to 
W-h jet production appeared 14411 
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partonic subprocesses in which the high px heavy boson 
is radiated from a more energetic quark (or a parton in 
general), as is illustrated in Fig. [TJ For the production 
of EW gauge bosons with mass Q ^ 100 GeV, the frag¬ 
mentation logarithm, \n{p\/Q^) ~ 4.6, when the boson’s 
transverse momentum pT ~ 1 TeV, and could be a po¬ 
tential source of large high order corrections. 




FIG. 1: Leading order QCD diagrams that lead to the frag¬ 
mentation logarithms in W and Z production. The diagram 
on the left (right) provides the fragmentation logarithm from 
quark (gluon) splitting. 


The leading contribution from the fragmentation log¬ 
arithm of the diagrams in Fig. [1] has the approximate 
form, 


cr^/a) 


27T 





dk^ 


f 

2 ^ W) 


( 1 ) 


where dj®, with parton flavor / = q,g, represents the 
lower order contribution to the production cross section 
without the quark to quark-gluon (or gluon to quark- 
antiquark) splitting in Fig. [U and Pf^q{z) is the LO 
splitting function for a parton of flavor / to fragment into 
a quark, e.g., Pg^q{z) = (l/ 2 )[ 2 ;^ + (l —z)^] with the color 
factor 1/2. To obtain the second line in Eq. ([T|), we use 
the fact that the fragmentation contribution is dominated 
by the large z region and the size of the LO splitting 
functions to a quark at large z is of 0{1). Since the 
factorization scale could be chosen from 0(Q) to 0{pt), 
we could estimate the size of the high order corrections 
from the LO fragmentation logarithms as 


0^1 f 

2tt \Qy 
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for pr ^ 1 TeV. That is, the higher order correc¬ 
tions from the fragmentation logarithms to the produc¬ 
tion of heavy EW gauge boson of mass Q ~ 100 GeV, 
[(ofs(p)/27r) ln(p|,/(5^)]™ with m > 1, should be under 
control perturbatively. It was pointed out in that 
the double logarithms from radiation of soft Z bosons 
can induce large higher-order corrections in the large jet 
Pt region oi Z+ jets production. This situation is dif¬ 
ferent from the case studied here where we are looking 


at inclusive Z boson production at large pt- The double 
logarithms are highly suppressed in our case. 

In this paper, we use a modified QCD factorization 
formalism for calculating the hadronic cross sections of 
W and Z production at large transverse momentum 
Pt- This formalism is similar to the one that we in¬ 
troduced for the low mass Drell-Yan cross section at 
large px in previous work (46l - l^ . In order to resum 
the a™ ln’"(p^/Q^)-type large logarithms, we re-organize 
the conventional fixed-order perturbative expansion of W 
and Z cross sections at large px into factored “direct” 
and “fragmentation” contributions, as demonstrated in 

Sec.in 

dCrAB^ViQ)X _ ^^AB^V{Q)X AB^V{Q)X 

dp^ dy dp^ dy dp^ dy 

All powers of the logarithmic ««ln(p|,/(5^) high order 
corrections are completely resummed into W and Z frag¬ 
mentation functions, and included into the “fragmen¬ 
tation” contributions, while the short-distance partonic 
hard parts for both “direct” and “fragmentation” con¬ 
tributions are free of the large logarithms, and can be 
systematically calculated order-by-order in powers of as ■ 
Using our modified factorization scheme, we calculate the 
next-to-leading order (NLO) predictions for W and Z bo¬ 
son production at high px at LHC energies, as well as at 
a future 100 TeV proton-proton collider. The perturba¬ 
tively calculated W and Z cross sections at large px show 
improved convergence, as well as moderate reduction of 
the scale variations compared to the conventional NLO 
perturbative expansion. 

The situation in W and Z production is very differ¬ 
ent from the fragmentation contributions to light hadron 
production, where the logarithms could run into the non- 
perturbative region, and also numerically different from 
the role of fragmentation contributions to the production 
of low mass Drell-Yan pairs or heavy quarkonia at high px 
[ttL [ 4 ^ . The key difference is the mass of the heavy EW 
gauge bosons. In order to produce a heavy EW gauge 
bosons of mass Q, as shown in Fig.[Tl the invariant mass 
of the fragmenting parton of momentum k should be suf¬ 
ficiently large, as > Q. Consequently, the radiation 
from the fragmenting parton is similar to the radiation 
from a heavy quark of mass Vk^ ^ Q, and is strongly 
suppressed for the phase space within the angle, Vl^/ko, 
of the fragmenting parton, the so-called “dead-cone” ef¬ 
fect [13 ■ In addition to the much smaller phase space for 
the radiation, the large virtuality of the fragmenting par- 
ton also sets the renormalization scale p for the strong 
coupling constant as{p) < as{Mz) 0.II8 with Z mass, 
Mz- With the current limit of collision energies, it is the 
combination of small as and the restricted phase space 
for the radiation that controls the size of the corrections 
from high order fragmentation logarithms. In this paper, 
we verify this conclusion by performing explicit all order 
resummation of the fragmentation logarithms. 

The rest of our paper is organized as follows. In 
Sec.in we introduce our modified factorization scheme 
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for EW gauge boson production at large pt and compare 
it with the conventional fixed-order perturbative expan¬ 
sion scheme. In Sec. lIIIl applying our modified factoriza¬ 
tion scheme, we present our predictions for the W and Z 
fragmentation functions, as well as our calculations of the 
LO and NLO cross sections of W and Z production at 
large pr at LHC energies and at a ^/s = 100 TeV future 
proton-proton collider. We also discuss the improvement 
of our modified factorization formalism over the conven¬ 
tional fixed-order perturbative expansion. Our summary 
and conclusions are presented in Sec. IlYl 


II. QCD FACTORIZATION OF VECTOR 
BOSON PRODUCTION AT LARGE pr 

The cross section for producing an on-shell EW gauge 
boson of momentum p^ can be factored systematically in 
QCD perturbation theory as [5lj | 

/ dx2fb{x2,p) 

J-Pert 

<^^ab^V{Q)X , „ ^ /Q^ 

X— -rr --.— {xi,x2,Q,PT,y,d-), (3) 

dprp cty 

under the usual assumption that the physically measured 
quantities Q and pp are both much larger than Aqcd- In 
Eq. ([3]), the variables Q, pr, and y are the mass, trans¬ 
verse momentum, and rapidity of the vector boson, re¬ 
spectively; and the symbol X stands for an inclusive sum 
over final states that recoil against the observed vector 
boson. The sum ^ runs over all parton flavors, and 
and are the corresponding PDFs, with the partons’ 
momentum fractions xi and X 2 , respectively; p repre¬ 
sents the renormalization and factorization scale, which 
is of the ord er of the energy exchange of the hard colli¬ 
sion: p ^ + p\. The function /dp'^dy 

in Eq. ([3]) is the short-distance hard part of the partonic 
scattering cross section. It can be calculated in perturba¬ 
tion theory in powers of the QCD coupling as{p). The 
first two terms of the perturbative expansion, dd‘^~^‘^ 
and including perturbative contributions up to 

O^aemC^s) and 0{aemCtl) respectively have been avail¬ 
able in the literature for sometime [§-[l0l|. The super¬ 
script “C” indicates the conventional fixed-order pertur¬ 
bative expansion. 

Beyond the leading order in when pr ^ Q, the 
perturbative functions d&^^J^y^Q^j^/dp'^dy in Eq. ([3]) 
can receive large high order corrections in powers of 
as hi(jj^/Q^) arising from the radiation of partons along 
the direction of the observed vector boson. Such large 
logarithmic corrections can be systematically resummed 
into the parton to vector-boson fragmentation functions, 
Df^y, as demonstrated for the case of the virtual pho¬ 
ton [^. The perturbative series for da^^J^yf^Q^^/dp'^dy 


can therefore be re-organized into two terms as in [47 


J-Pert 

^^ab^V{Q)X , „ , 

- -T^ - {Xi,X2,Q,PT,y,d') = 

clprp cLy 


^^ab^V(Q)X 

dpi, dy 


{xi,x2,Q,PT,y;p,PD) 


\Jrj^ 

T- Frag 

^^ah^V{Q)X 

dp\ dy 


{xi,X2,Q,PT,y,P,PD) , (4) 


where the superscripts “Dir” and “Frag” represent the 
“direct” and the “fragmentation” contribution, respec¬ 
tively. The latter includes the perturbative fragmenta¬ 
tion logarithms, and it can be further factored as [fi^ . 


^^ab^V{Q)X 

dpidy 


E 


dz 

72 


dda 


b^cX 


dpl,^ dy 

2 , r\2 


{xi,x2,Pc;yD) 


P£),Q ), 


(5) 


where p>^ = p^ jz, with defined to be at = 0, 
corresponding to the approximation <C made for 
ddab^cX/dp^^ dy in Eq. ([5]), which is a short-distance 
hard part for partons of flavors a and b to produce a par- 
ton of flavor c and momentum p^. The fragmentation 
scale pd ^ pt is introduced to separate the direct and 
fragmentation contributions in Eq. o, and its depen¬ 
dence would be cancelled if both contributions include 
all order corrections [47l |. 

The fragmentation contribution in Eq. m shares the 
typical two-stage generic pattern of the fragmentation 
production of a single particle at large transverse mo¬ 
mentum Pt (much larger than its mass): the production 
of an on-shell parton of flavor c at the distance scale 1/pt, 
convoluted with a fragmentation function that includes 
the leading logarithmic contributions from the “running” 
of the distance scale from l/pn ^ ^/pt to l/poo ~ 1/Q- 
When Q 7^ Aqcd , which is the case for the W and Z bo¬ 
son, the parton to vector-boson fragmentation function 
Df^v is perturbative, and so is the whole resummed 
fragmentation contribution, in Eq. ® M- 

The direct contribution in Eg. (l4l) is perturbatively cal¬ 
culable in a power series of as (47|, and it is defined as. 


°‘^ab^V(Q)X 

dp\ dy 


pdPert 

'^^ab^V{Q)X 

dp?p dy 


da 


F—Asym 
ab^V{Q)X 

dp^ dy ’ 


( 6 ) 


where is the conventional fixed-order perturbative 

QCD calculation. Moreover, with the super¬ 

script “Asym” referring to an “asymptotic” contribution, 
is simply the perturbative expansion of the resumed frag¬ 
mentation contribution, in Eq. ([5]) , to the same or¬ 

der in powers of as as dd^®''* in Eq. (jS]). It is effectively 
a subtraction term to systematically remove from dd^®"'* 
all fragmentation logarithms which have been resummed 
into the fragmentation contribution, dd^’’®® in Eq. ©■ 
The subtraction avoid double counting order-by-order in 
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powers of as- Consequently, the direct contribution is 
free of the large fragmentation logarithms, while it still 
keeps all non-logarithmic terms from the physics between 
the scales px and Q. Note that since 

the fragmentation contributions start at 0{aemCt^)- 
By substituting Eqs. o and ([5]) into Eq. ([3]), we ob¬ 
tain our modified factorization formalism for heavy EW 
vector boson production in hadronic collisions. 


dcr ab^v{Q)x 
dp^ dy 


^^AB^V(Q)X 

dp^ dy 


da 


Frag 

AB^V(Q)X 
dp^ dy 


(7) 


where da^^’^ is factored in the same way as da in Eq. © 
with replaced by d(T^“' defined in Eq. ®, and 

is 


da 


Frag 

AB^V(Q)X 

dp^ dy 



where the fragmentation functions Dc^v{zt re¬ 
sum all fragmentation logarithms, and are par- 

tonic hart parts for producing an on-shell parton of flavor 
“c” and momentum p^, which are independent of the spe¬ 
cific vector boson produced. Actually, they are the same 
as the perturbative coefficient functions for producing a 
light hadron, such as pion, and are available for both the 
LO and NLO in powers of as in the literature [sS j. 

The separation between the direct and the fragmen¬ 
tation contribution in Eq. ([7|) depends on the definition 
of the parton to vector-boson fragmentation functions. 
Different definitions of the fragmentation functions cor¬ 
respond to a different scheme to split the conventional 
fixed-order perturbative expansion into the “direct” plus 
“fragmentation” contributions. A scheme choice for the 
fragmentation functions should also fix the asymptotic 
contribution and, correspondingly, the direct contribu¬ 
tion. The sum of these two terms in our modified factor¬ 
ization formalism in Eq. © should not be very sensitive 
to the scheme choice. 

The key difference between our modified factorization 
formula in Eq. © and the conventional factorization for¬ 
mula in Eq. Q resides in the way the large logarithmic 
contributions from final-state parton splitting are han¬ 
dled. Instead of one perturbative series in powers of as 
in the conventional approach, we have two perturbative 
expansions in our modified factorization formula: one for 
the direct and one for the fragmentation contribution. 
All coefficient functions in the new perturbative expan¬ 
sions are free of large logarithms. The large perturba¬ 
tive logarithms in the conventional expansion in powers 
of as are systematically resummed into the fragmenta¬ 
tion functions. As a result, the perturbative expansion 
of the partonic hard parts in our modified factorization 


approach has better convergence properties than the con¬ 
ventional fixed-order expansion. In addition, because of 
the reorganization and resummation of the logarithms, 
both the LO and NLO cross sections in our new approach 
include a tower of logarithmically enhanced high order 
corrections from the conventional fixed-order perturba¬ 
tive expansion. Eor example, for the LO cross section, 
with the lowest order partonic hard parts retained, the 
difference between our modified approach in Eq. © and 
the conventional fixed-order approach in Eq. is 



da 


D-LO 

ab^V{Q)X 


dp^ dy 


da 


C-LO 

ab^V{Q)X 


dp^ dy 


E 5 


dz da 


F-LO 
ab—^cX 

dpl^ dy 

<'dz dal 


D. 


c^V{Q) 


Ct 

F-LO 

ab—>^cX 


a.b.c 




D 


c-*-V(Q)i 


(9) 


where the LO relation, = dd^ is used. 

Similarly, at NLO, the difference is 


E JdxifaJdx2fi^ 

a,b 


t-D-NLO 

^"ab-pViQ)X 


.-C-NLO 


dpy dy dp^ dy 


E 


^"ab^cX jj 

z^ dpl^ dy 


c^ViQ) 


E^ 

a,b,c 


dxifi 


dx 2 fn^ 


dz 
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dd 


F-LO 

ab—^cX 


dpl^dy ^c^viQ)\ 


-D 


(0) 


-k 


dd 


CT 

F-NLO 

ab—^cX 


dd 


F-LO 

ab—^cX 


dpl^ dy dpl^ dy 


D 


c^V(Q) 


• ( 10 ) 


In deriving Eq. m, we use the definition in Eq. ® and 
the following relation. 


F-Asym-NLO 

^^ab^ViQ)X 

dp^ dy 


^ ddlb^?x ^( 0 ) 
^ J z^ dpl^ dy 


( 11 ) 


where is the fragmentation function for a sin¬ 

gle parton of flavor c to fragment into a vector boson of 
invariant mass Q, expanded perturbatively to 0{ae.mOcl)- 
In comparison with a conventional fixed-order calcula¬ 
tion at the same orders, Eqs. ® and (fTUl) show that our 
modified factorization formalism in Eq. © includes the 
all orders resummation of final-state fragmentation loga¬ 
rithms, plus additional short-distance high order contri¬ 
butions. The first term on the right-hand side of Eq. m, 
comes from the resummation of leading logarithms to all 
orders, and the second term is from the 0{a^) correc¬ 
tions to the partonic hard parts. To the order in which 
we work, the full correction is the sum of the two. In 
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our presentation of numerical results in Sec. nil the sum 
of the two terms in Eq. (fTUl) is referred to as the type-B 
correction to the NLO term of the conventional fixed- 
order perturbative calculation, while the first term alone 
is denoted the type-A correction. 

The resummation of the fragmentation logarithms into 
the single parton to a vector-boson fragmentation func¬ 
tions is achieved by solving the inhomogeneous evolution 
equations (s^ . 

+ (^) X! / —P<^^d{^,a8)Dd^v{z',n\,\Q‘^), 

d 

( 12 ) 


where c, d = g, g, g. The ambiguity in defining the frag¬ 
mentation function is connected to the renormalization 
of its perturbative ultraviolet (UV) divergence, and the 
choice of the fragmentation scale, go- In Eq. (HI, 
the evolution kernels Pc-s-d are infrared (IR) safe, eval¬ 
uated at a single hard scale, go, and calculated pertur- 
batively as a power series in a^. These kernels are in¬ 
sensitive to the specific vector meson produced, and are 
the same as the evolution kernels for pion or other light 
hadron fragmentation functions. They are the same as 
the splitting functions of the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution equations [sj at the 
leading order. 

The inhomogeneous term in the evolution equations 
can also be calculated perturbatively, and in general it 
has power correction terms of the form owing 

to the mass of the vector boson, and therefore, it de¬ 
pends on the nature of the vector boson produced. In 
the invariant mass cutoff scheme [s^, , the lowest or¬ 
der quark-to-virtual photon evolution kernel was 

derived in [^ . With a simple generalization, we com¬ 
pute the inhomogeneous evolution kernel for a quark to 
fragment into an EW gauge boson of invariant mass Q, 
obtain, 



{\9lr 




l + il-zf 





), (13) 


vanishes at the lowest order 

-f^%y{z,e-,Q^) = 0, (14) 

because the gluon does not interact directly with the EW 
gauge bosons. 

The choice of factorization scheme is not unique. We 
use the invariant mass cutoff scheme [s^, 1^ . Some 
choices, such as the modified minimum subtraction (MS) 
scheme, may not respect the mass threshold when 
and lead to negative fragmentation functions j52l 
QCD corrections to the lowest order parton to vector- 
boson splitting function can be evaluated in prin¬ 

ciple order-by-order in as ■ 

If Q Aqcd, the parton to vector-boson fragmenta¬ 
tion functions are completely perturbative, as is the case 
for our study here. The lowest order parton to vector- 
boson fragmentation functions are 





i\9lr 




/ Ck-em \ 

V^j 


D%{z,g^j,-Q^) = Q. 


\Q^ J 


ZfJ-lJ. 


,(15) 

(16) 


We must specify a boundary condition in order to solve 
the evolution equations in Eq. (ini). A natural boundary 
condition following the mass threshold constraint is |4^ . 


D,^y{z,gl<Qyz-Q^) = 0, (17) 

for any flavor c, if we choose the invariant mass cutoff 
scheme for the fragmentation functions. 

To conclude this section, we emphasize that our modi¬ 
fied factorization formalism in Eq. o effectively reorga¬ 
nizes the single perturbative expansion of conventional 
QCD factorization, in Eq. ©, into two perturbative ex¬ 
pansions in Eq. o, plus the perturbatively calculated 
evolution kernels for the fragmentation functions. The 
main advantage of this reorganization is that the new per¬ 
turbative expansions are evaluated at a single hard scale 
and are free of large logarithms. As shown in Sec. IIIII 
the ratios of the NLO over the LO contributions in the 
new perturbative expansions are smaller than the ratios 
evaluated in the conventional approach. 


where is the invariant mass of the parent quark and 
is identified as and the 0-function is a consequence 
of the mass threshold. Here, g'^^ are the EW cou¬ 
plings between quarks and the EW gauge bosons, with 
g^'^} = {l/(-\/2siu), 0} for the W boson, and 
{9l‘^, 9r^} = {(1/2 - 2s^/3)/(su,Cu,), -2su,/(3cu,)} and 
{(-1/2 -I- s^/3)/(su,Cu,), +Syj/{‘icyj)} for the Z boson to 
couple to the up and down-type quarks respectively. We 
us Sw and to represent the sine and cosine of the weak- 
mixing angle. The gluon to vector-boson evolution kernel 


III. NUMERICAL RESULTS 

We present our numerical results in this section. We 
choose the Gp parametrization scheme for the EW 
couplings with Mw = 80.385 GeV, Mz = 91.1876 GeV, 
Mt = 173GeV, and Gp = 1.166379 x IQ-'^GeV^ [s^. 
We assume a diagonal Cabibbo-Kobayashi-Maskawa ma¬ 
trix of the SM in the calculation for simplicity. We use 
METAl.O PDE [ 5 ^ of the proton, which is a NNLO 
PDF set that combines information from the CTIO [s^ . 
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MSTW2008 [H, and NNPDF2.3 [HI PDF sets, and 
three-loop running of the QCD coupling constants with 
Nf = 5 active quark flavors and as{Mz) = 0.118. 

A. Fragmentation functions 

At LO only quarks can fragment into EW vector 
bosons, for which the fragmentation functions are pro¬ 
portional to the square of the EW couplings. As seen 
in Eq. (HSl), for example, the up-type quark can frag¬ 
ment into a W~^ boson at LO, but not into a W~ boson. 
The fragmentation functions are the same for quark and 
anti-quark to virtual-photon or Z boson. Similarly frag¬ 
mentation functions are equal for up-type quarks to W'^ 
and down-type quarks to 1F“. The gluon can fragment 
into W or Z bosons only at higher orders through in¬ 
termediate quarks, thus the fragmentation functions are 
suppressed. 

In Fig. [21 we show the momentum fraction 2 depen¬ 
dence of the LO and resummed parton to vector-boson 
fragmentation functions for the Z and W~^ bosons. In 
each figure we also plot the parton to virtual-photon frag¬ 
mentation functions with the same mass. We choose a 
typical fragmentation scale of 800 GeV and show results 
for the u quark as an example. The fragmentation func¬ 
tions have a peak in 2 ;, a consequence of the vector bo¬ 
son’s mass threshold. The height of the peak is propor¬ 
tional to the square of the corresponding EW couplings. 

In Fig. 121 we show the resummed corrections to the 
LO fragmentation functions in percentage. As expected, 
the resummed contributions increase the fragmentation 
function at small and moderate z values while they re¬ 
duce it for large z values. The observed differences of 
the fractional corrections for different vector bosons are 
caused by the difference of their EW couplings through 
the mixing of singlet and non-singlet evolutions. We fur¬ 
ther show the evolution of the fragmentation function 
for a u quark into a Z boson with momentum fraction 
z = 0.04 in Fig. [H The scale dependence of the LO 
fragmentation functions are dominated by the logarith¬ 
mic term in Eq. (USD. The resummed corrections here are 
always positive. 


B. W and Z boson production at large pr in the 
conventional expansion scheme 

In this subsection we present the cross sections of W 
and Z boson production at large pr at LO and NLO 
accuracy in the conventional fixed-order perturbative ex- 
pansioUjdefined in Eq. ®, based on the analytical re¬ 
sults in [§-[1(TI|. We plot in Fig.jSlthe pr dependence of the 
double-differential cross sections of the vector bosons at 
a central rapidity, y = 0- We show the results at y/s = 8 
and 14 TeV at the LHC, and at y/s = 100 TeV at a fu¬ 
ture proton-proton collider, with increasing range of px 
as allowed by the higher collision energy. We focus on the 




FIG. 2: The z dependence of the LO and resummed fragmen¬ 
tation functions of a u quark into a Z and a boson for 
a fragmentation scale yn ~ 800 GeV. For reference, we also 
plot the fragmentation function of a virtual photon with the 
same mass. 


large px region (> 500 GeV) where the perturbative log¬ 
arithms, oc \n{px/Q), are large. We set the factorization 
and renormalization scales, as well as the fragmentation 
scale, equal to the transverse mass mx = y/Q^ + Px^ 
unless specified otherwise. The scale dependence of our 
numerical results is discussed later. As shown in Fig. |5| 
the cross sections decrease quickly with increasing px at 
the LHC, caused by the fast drop of the partonic flux at 
large partonic momentum fraction x. With a larger phase 
space available at the 100 TeV collider, the px spectrum 
decreases more slowly. The NLO QCD corrections can 
be as large as 50% of the LO cross sections for the LHC 
energies, and reach 60% for the 100 TeV collider. 

To quantify the size of contributions from the large 
fragmentation logarithms we compare the full NLO cor¬ 
rections with the part of the corrections proportional to 
the fragmentation logarithms, equal to the fragmentation 
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FIG. 3: Percentage change of the resummed to the LO frag¬ 
mentation functions for a u quark to (a) Z and (b) W'^ boson 
for a fragmentation scale hd ~ 800 GeV. 


contribution in Eq. (jS]) expanded to the same order in as- 
In Fig. ini we plot the normalized full NLO corrections, 

(solid lines), along with the 
normalized contributions proportional to the fragmenta¬ 
tion logarithms (dashed lines). The full NLO corrections 
are large and increase quickly for vector-boson produc¬ 
tion at the LHC as px approaches its maximum. The 
behavior differs slightly for different vector bosons owing 
to the different masses and different contributing PDFs. 
The fragmentation pieces contribute about one third of 
the full NLO corrections for the px values considered. 
With the larger phase space available for radiation, the 
NLO fragmentation contributions are more significant at 
the 100 TeV collider. Their magnitude is as large as 30% 
of the LO cross sections or about one half of the full NLO 
corrections. The fragmentation contributions are almost 
insensitive to the value of px for the 100 TeV collider, 
and they stay pretty much constant when pt ^ 1 TeV 



liD [GeV] 

FIG. 4: Evolution of the LO and resummed fragmentation 
functions for a u quark into a Z boson with momentum frac¬ 
tion 2 = 0.04. 


for both W and Z bosons. This result comes about be¬ 
cause the logarithmic enhancements from \ti(j)x/Q) are 
also modulated by the 2 dependence of the hard partonic 
cross sections in the convolution with the fragmentation 
functions. 

Equation (jS]) expresses the fragmentation contribution 
as a convolution of two PDFs and one fragmentation 
function, and the hard scale of the partonic scattering 
is proportional to the combination xi x^jz. The factored 
fragmentation contribution should likely be dominated 
by the kinematic region where xi and X 2 are small, while 
2 is relatively large, owing to the competition of two 
steeply falling PDFs against one fragmentation function 
when the momentum fractions, a; or 2 increase. To better 
understand the dominant region of partonic scattering, 
we introduce a cutoff, 2cut < 1, to limit the 2 -integration 
in Eq. ([5]) to 2niin < Z < 2cut, instead of 2niin ^ z < 1. 
In Fig. [71 we plot fQj. ^ boson production as a 

function of the cutoff, 2cut, normalized by the full contri¬ 
bution with 2cut = 1, to show the fractional contributions 
from a limited phase space in z. As expected, the contri¬ 
butions at large 2 dominates in the fragmentation pieces. 
For example, the domain 2 > 0.5 contributes over 90% of 
the fragmentation cross sections for all collision energies 
shown in Fig. [71 The effective range of 2 is even more 
limited to the larger 2 -values for production at the LHC 
energies owing to the smaller phase space for a given px- 

In Fig. El we plot the LO and NLO contributions to the 
double-differential cross sections of Z boson production 
at j/ = 0, evaluated in the conventional fixed-order per¬ 
turbative expansion and normalized by the conventional 
LO cross sections. The uncertainty bands in Fig. El are 
calculated by varying the renormalization and factoriza¬ 
tion scales, pr = P/ = M, upward and downward by a 
factor of two, as an estimation of the remaining higher- 
order contributions. There are moderate reductions of 
the scale variation from LO to NLO for the LHC energy, 

















FIG. 5: The LO and NLO predictions for the double¬ 
differential cross sections of Z and boson production at 
large transverse momentum, pT, in the conventional fixed- 
order perturbative expansion, at -y/s = 8 and 14 TeV at the 
LHC, and 100 TeV at a future proton-proton collider. The 
range of pr is determined by the available phase space of the 
collisions. 

while the LO result at the 100 TeV collider shows a scale 
variation that is too small due to accidental cancellations 
between the renormalization and the factorization scale 
dependence. One may also notice that the NLO bands 
do not overlap the LO bands, suggestive that the scale 
variations for the conventional fixed-order perturbative 
expansion may underestimate the true theoretical uncer¬ 
tainties from the remaining high order corrections. 

C. W and Z boson production at large pt in the 
modified expansion scheme 

In our modified factorization scheme, we split the con¬ 
ventional fixed-order perturbative series into the “direct” 
and “fragmentation” contributions, as shown in Eq. 0, 



FIG. 6: The full NLO corrections in the conventional fixed- 
order expansion scheme, and the corresponding contributions 
proportional to the fragmentation logarithms, both normal¬ 
ized to the conventional LO cross sections, are plotted as a 
function of pr for Z and boson production in proton- 
proton collisions at = 8, 14 and 100 TeV. 


and we calculate them to NLO in this paper. In Fig.|9l we 
plot both the LO and NLO “direct” contributions to the 
double-differential cross sections of Z boson production. 
The results are presented as a ratio , normalized by the 
conventional LO cross sections. Note that as mentioned 
earlier the “direct” contribution at LO equals the con¬ 
ventional LO cross section. Similarly, in Fig. [TOl we plot 
the corresponding LO and NLO “fragmentation” contri¬ 
butions, normalized by the same conventional LO cross 
sections. In our modified factorization scheme, there is 
a fragmentation scale po, in addition to the renormal¬ 
ization scale p,r and factorization scale p /, which are the 
same as those in the conventional fixed-order factoriza¬ 
tion scheme. Different choices oi fijy effectively move 
some finite perturbative contributions between the “di¬ 
rect” and the “fragmentation” contributions in Eq. ©■ 
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■^cut 



FIG. 7: Fractional contribution of for Z boson pro¬ 

duction as a function of the cutoff, Zcut, of the upper limit 
of a-integration in Eq. ((8ll, at pr = 500 and 1000 GeV and 
y = 0. 


We set the default choice for all three scales equal to uit- 
The uncertainty bands in Figs. [9] and [TO] are calculated 
by varying the renormalization and factorization scales, 
Mr = M/ = Mi upward and downward by a factor of two, 
as an estimation of the remaining higher-order contribu¬ 
tions. The effects from varying the fragmentation scale, 
flu, are discussed later. 

Comparing the relative size of NLO corrections in our 
modified factorization scheme with those of the conven¬ 
tional fixed-order expansion, we see in Figs. [51111 and [TUI 
that the NLO K-factors in the conventional expansion 
are significantly larger than those evaluated in the mod¬ 
ified scheme. The improved convergence of our modified 
expansion scheme indicates a better control of unknown 
higher-order corrections. In Fig. 1101 we observe signifi¬ 
cant reduction of the scale variation for the fragmentation 
contributions at NLO as a consequence of the inclusion 
of the NLO partonic hard parts. 



Pt,z [GeV] 



Pt,z [GeV] 

FIG. 8: Conventional LO and NLO predictions on double¬ 
differential cross sections of Z boson production at 14 TeV 
and 100 TeV, including the central predictions and the scale 
variations, all normalized to the central predictions of the 
conventional LO results. 


As expressed in Eq. ©, the full NLO predictions in 
our modified factorization scheme are equal to the sum 
of the NLO “direct” and “fragmentation” contributions. 
They differ from the conventional NLO results by ad¬ 
ditional higher-order contributions, or the type-B correc¬ 
tions, defined in Eq. (flUl) . In Fig.[TTl we plot the full NLO 
predictions in our modified factorization scheme for the 
double-differential cross sections of Z boson production, 
together with the conventional NLO results, where we 
show both the central predictions and the scale variation. 
Both are normalized by the central values of the conven¬ 
tional LO cross sections. Note that the scale-variation 
(blue) band of the modified factorization scheme lies al¬ 
most entirely within the (red) band of the conventional 
scheme. At NLO accuracy, the predictions from the con¬ 
ventional fixed-order perturbative expansion show about 
20 % uncertainty from the scale variation, which increases 
withpp, while as expected, the full NLO predictions from 
our modified factorization scheme show a reduction of the 
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Pt,z [GeV] 



Pt,z [GeV] 

FIG. 9: LO and NLO predictions for the “direct” contribu¬ 
tions to double-differential cross sections of Z boson produc¬ 
tion in proton-proton collisions at s/s = 14 TeV and 100 TeV, 
respectively, with uncertainty bands from the scale variations. 
The predictions are normalized by the central values of the 
conventional LO cross sections. 



Pt,z [GeV] 



Pt,z [GeV] 

FIG. 10: LO and NLO predictions for the “fragmentation” 
contributions to double-differential cross sections of Z boson 
production in proton-proton collisions at = 14 TeV and 
100 TeV, respectively, with uncertainty bands from the scale 
variations. The predictions are normalized by the central val¬ 
ues of the conventional LO cross sections. 


uncertainty, because of the inclusion of the type-B correc¬ 
tions. The reduction is stronger for the NLO predictions 
at y/s = 100 TeV, where the fragmentation contributions 
are significantly larger, as shown in Fig. [6l 

The central values of the NLO predictions from our 
modified factorization scheme are larger than the NLO 
predictions from the conventional factorization scheme, 
as shown in Fig. [TTJ This result is an immediate conse¬ 
quence of the type-B corrections in Eq. (fTOll . To fur¬ 
ther quantify this difference, we plot the type-B cor¬ 
rections for both the Z and boson production in 
Fig. nil along with the type-A corrections, as speci¬ 
fied in the discussion immediately after Eq. (fTTll . Both 
type-A and type-B corrections in Fig. [T^j are shown as 
ratios to the conventional LO cross sections, Rcor. = 
(^D-NLO ^F-NLO _ ^C-NLO)/^C-LO ^ TT ulih e the type- 

B corrections, which are positive and defined to be the 
total corrections at NLO, the type-A corrections are neg¬ 


ative and are of a few percent of the LO cross sections. 
As shown in Eq. the type-A corrections are propor¬ 
tional to the difference between the fully-resuinmed and 
the LO perturbative fragmentation functions. The fully 
resummed fragmentation functions have a softer z depen¬ 
dence when compared to the LO piece of the functions, 
and naturally, the difference is negative for the large z re¬ 
gion, while it becomes positive and larger as z decreases, 
as shown in Fig. |31 The negative type-A corrections, as 
shown in Fig. [T^l arise from the fact that the fragmenta¬ 
tion contributions are dominated by the large z region, 
as discussed above. Numerically, the type-A corrections 
are within —6% for all cases shown in Fi g. 1121 similar to 
the results for virtual photon production . At a given 
Pt, the absolute size of the corrections is larger for pro¬ 
duction at a lower owing to the smaller phase space 
for radiation and larger effective z values. The net type- 
B corrections are positive, despite the negative type-A 
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Pt,z [GeV] 



FIG. 11: The full NLO predictions in the modified factoriza¬ 
tion scheme for double-differential cross sections of Z boson 
production at 14 TeV and 100 TeV, including the central pre¬ 
dictions and the scale variations, together with the conven¬ 
tional NLO results, all normalized to the central predictions 
of the conventional LO results. 


corrections, because the fragmentation corrections from 
NLO hard matrix elements are positive, bringing the to¬ 
tal corrections up by 6-10%. This outcome is understand¬ 
able since the NLO corrections to parton or jet produc¬ 
tion at hadron colliders, which share the similar hard 
matrix elements, are usually positive and large. 

As shown in Fig. [T^l the final numerical differences 
between the NLO predictions, calculated in our modi¬ 
fied factorization scheme and the conventional fixed-order 
perturbative expansion are about 5% at pt = 500 GeV, 
and increase or decrease slightly for the case of the LHC 
or a 100 TeV proton-proton collider. The corrections 
are similar for Z and W'^ boson production. Based on 
the relatively weak px dependence of the corrections in 
Fig. [T21 we expect that the NLO predictions from our 
modified factorization scheme will not differ appreciably 
for the overall shape of the large px spectrum of the W 



FIG. 12: Differences between the full NLO predictions in the 
modified factorization scheme and the conventional NLO pre¬ 
dictions, for double-differential cross sections of Z and 
boson production at 8, 14, and 100 TeV, in each sub-figure 
with increasing range of px- 


or Z at the LHC derived from the conventional NLO cal¬ 
culations. At the 100 TeV proton-proton collider, our 
NLO predictions for the large px spectrum of both the 
W and Z boson production could be slightly softer than 
the NLO spectrum evaluated from conventional fixed- 
order perturbative calculations without resummation of 
the fragmentation logarithms. It is the resummation of 
the large fragmentation logarithms that makes our mod¬ 
ified factorization formalism more stable and converge 
faster perturbatively. 

To examine the fragmentation scale dependence of the 
NLO predictions of our modified factorization scheme, we 
increased and decreased the scale by a factor of two from 
its default value, po = rriT, while fixing the renormal¬ 
ization and factorization scales. Our NLO results in the 
modified factorization scheme show a very weak depen¬ 
dence on px), at a level of 1%. By comparing the frag¬ 
mentation scale dependence between the “direct” and the 
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“fragmentation” contribution to the production cross sec¬ 
tions, and in Eq. o, we observe that the leading 
dependence on in cancels that in by definition. 
The sub-leading dependence on fio from QCD evolution 
cancels within between the resummed corrections to 
the fragmentation functions and the NLO corrections to 
hard matrix elements. Thus the remaining dependence 
on fj,D in the sum of q-D-nlo ^f-nlo jg negligible. 

In conclusion, within the modified factorization scheme 
we split vector-boson production at large pt into two sep¬ 
arate perturbation series, both of which show better con¬ 
vergence compared to the conventional fixed-order calcu¬ 
lations. By including some higher-order corrections (be¬ 
yond NLO in conventional fixed order) for the fragmen¬ 
tation contributions we reduce the scale variations and 
obtain a perturbative expansion with better control of 
the theoretical uncertainties, especially for vector-boson 
production at the 100 TeV collider. 


D. Discussion of the Numerical Results 

In Table |T] we summarize our total cross sections for 
vector-boson production, integrated over the full rapid¬ 
ity range, with pry greater than 800 GeV. We provide 
predictions based on conventional LO and NLO in the 
second and third columns, and the NLO predictions us¬ 
ing the modified factorization scheme in the last column, 
for the LHC at 13 TeV and for a future proton-proton 
collider at 100 TeV. We choose a pr threshold of 800 GeV 
for illustration. As shown in Table HI the total cross sec¬ 
tion for Z boson production at y/s = 13 TeV can reach 
120 fb with a scale uncertainty about ±8%. The scale 
uncertainty is about ±5% for Z boson production at the 
100 TeV collider. Uncertainties due to variations of the 
fragmentation scale po are not included since they are 
much smaller. These residual scale variations should be 
further reduced either by soft gluon resummation Hi- 
1^ or ongoing NNLO calculations. Other uncertainties 
include PDF uncertainties of about ±4% for Z or W bo¬ 
son production at 13 TeV, and about ±1% for the 100 
TeV collider, obtained using the META PDF prescrip¬ 
tion proposed in The PDF uncertainties are much 
larger at 13 TeV because the large-x region of the PDFs 
is being probed, a situation that could be improved once 
more LHG data are included in the global analysis of 
PDFs. For vector-boson production at large pr at the 
LHG, the electroweak corrections due to Sudakov loga¬ 
rithms are also significant [lll - [l8j| and must be included 
when comparison is made with data. 

Measurements exist at y/s = 7 TeV at the LHG of 
the Z boson px spectrum in the leptonic channel @ and 
of the Z or W boson cross sections above a certain pr 
threshold in the hadronic channel Q- The highest effec¬ 
tive Pt of the vector boson is limited to about 300 GeV. 
At 13 TeV, according to Table HI ~2400 Z bosons are 
expected with px above 800 GeV in the di-muon and di¬ 
electron channel, assuming an integrated luminosity of 


CI'pt> 800 GeV 

VB 

C-LO 

C-NLO 

NLO [modified] 


Z 

74.1 

117.41)?;° 

120.5l?o^4 

13 TeV [fb] 

w+ 

126.2 

199.41?°;) 

204.41??;) 


W- 

55.8 

90.2l°;« 

92.71);) 


z 

11.48 

19.681?;“ 

20.161);“ 

100 TeV [pb] 

w+ 

15.08 

26.231?;?) 

26.861?;)) 


W- 

10.50 

18.181?;)) 

18.611);“ 


TABLE I: Cross sections integrated over y and px > 800 GeV 
with scale variations, for W and Z boson production at 13 
TeV (in fb) and at 100 TeV (in pb). The scale uncertain¬ 
ties are calculated by varying pr = Pf = p hy a, factor of 
two in both directions around mx- The uncertainties due to 
variations of po are not included since they are much smaller. 

300 fb~^. The statistical uncertainties are estimated to 
be smaller than the scale variations in Table HI Thus pre¬ 
cise tests of various QCD predictions require good control 
of the experimental systematical uncertainties at 13 TeV 
and 100 TeV. The situation is similar for W boson pro¬ 
duction where we gain a larger branching ratio in the 
leptonic channel but suffer from lower efficiency in the 
event reconstruction. 


IV. SUMMARY 

We introduced a modified factorization scheme for 
evaluating the transverse momentum spectrum of the SM 
W and Z boson production at hadron colliders. In this 
new scheme we re-organize the conventional fixed-order 
QCD perturbative expansion into two separate pertur¬ 
bative expansions, corresponding to “direct” and “frag¬ 
mentation” contributions. The fragmentation piece has 
a typical two-stage pattern, production of an on-shell 
parton convoluted with the perturbative fragmentation 
functions of the W and Z boson, as illustrated in Fig. [T] 
When p\ the large perturbative fragmentation 

logarithms of Iniji^/Q^) are resummed into the frag¬ 
mentation functions by solving their evolution equations. 
Consequently, the short-distance partonic hard parts for 
both the “direct” and the “fragmentation” contributions 
are free of these large fragmentation logarithms, and the 
overall convergence of the re-organized perturbative ex¬ 
pansions is improved, in particular, for the region where 

In our explicit NLO calculations, the fragmentation 
logarithms make up a large portion of the conventional 
NLO K-factors at high px at the LHC and especially 
at a future 100 TeV collider. Our improved NLO 
predictions retain all ingredients in the conventional 
NLO calculations (up to 0{aemO:l)), but they include 
partial higher-order corrections. In comparison with 
the conventional fixed-order perturbative expansion, the 
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modified NLO predictions show a moderate reduction 
of the scale variation at large PT^ and improved conver¬ 
gence. The improved NLO predictions are about 5% 
higher in the normalization of the large-pT spectrum, 
but they provide only modest changes in the shape 
of the spectrum. In subsequent research, it would be 
desirable to extend the analysis to NNLO. Our modified 
factorization scheme could be applied to the production 
of other heavy particles, such as the Higgs boson and 
top quark, at large pr when the px of the produced 
state is much larger than its mass. 
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